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Abstract. We report on the detection of the lowest CO J = 1 — transition in the powerful high-redshift radio 
galaxy 4C 60.07 at z = 3.79. The CO emission is distributed in two spatially and kinematically distinct components 
as was previously known from the observations of the higher excitation CO J = 4 — 3 line. The total molecular 
gas mass in 4C 60.07 inferred from the CO J = 1 — emission is M(H2) ~ 1.3 X 10^^ Mq , sufficient to fuel the 
inferred star- formation rate of ~ 1600 Mq yr~^ for 10* yrs. From our high-resolution CO J = 1 — VLA maps we 
find the dynamical mass of 4C 60.07 to be comparable to that of a giant elliptical at the present time. A significant 
fraction of the mass is in the form of molecular gas suggesting that 4C 60.07 is in an early state of its evolution. 
The merging nature of 4C 60.07 along with its large dynamical mass imply that this system is a giant elliptical 
caught in its formative stages. 
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1. Introduction 

High-redshift radio galaxies (HzRGs) are amongst the 
most luminous objects known, and are believed to serve as 
tracers of the peaks of the primordial density field around 
which giant elliptical galaxies and clusters of galaxies form 
(Kauffmann 1996; West et al. 1994). In the radio, HzRGs 
typically display a double-lobe morphology and large radio 
luminosities {PnsMHz ^ 5x 10'^^ ergs^^ Hz^^), indicat- 
ing a highly active black hole. 

Recently, evidence has been mounting that HzRGs are 
massive starburst galaxies. This has come about from sub- 
millimetre detections of a number of HzRGs, implying 
large rest- frame far-IR luminosities [LpiR — IQ^^ Lq ) 
powered by intense star formation [SFR ~ 1000 — 
2000 Mq yr^i - Dunlop et al. 1994; Hughes et al. 1997; 
Ivison et al. 1998; Archibald et al. 2001). In a recent 
SCUBA survey of seven HzRGs and their surroundings, 
Stevens et al. (2003) not only found the star forma- 
tion in the radio galaxies themselves to be extended on 
several tens of kilo-parsec scales but also found one or 
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more previously undetected submm sources in the vicinity 
(50 — 250 kpc) of more than half of the targeted objects. It 
is difficult to sec how the Active Galactic Nucleus (AGN) 
could power the far-IR luminosity on ^ 10 kpc scales, and 
a massive starburst seems to be the natural explanation. 
Indeed, adequate "fuel" for such large star formation rates 
has been found in the four HzRGs which have been de- 
tected in CO to date (Papadopoulos et al. 2000; De Breuck 
et al. 2003a; De Breuck et al. 2003b). These observa- 
tions revealed the presence of massive (^ 10"'^^ M© ) reser- 
voirs of molecular gas, enough to fuel a ~ 1000 M© yr~^ 
starburst for ~ 10^ yr, and in in half of the cases the 
CO emission was found to be extended on tens of kpc 
scales (Papadopoulos et al. 2000; De Breuck et al. 2003a). 
Similar large molecular gas masses distributed in clumps 
on tens of kilo-parsec scales has been found in a number 
of QSOs at high redshifts (Carilfi et al. 2002a; CariUi et 
al. 2002b). In general, HzRGs have the advantage over 
quasars that they are not gravitationally lensed since 
they are usually selected on the basis of extended lobe- 
emission whereas quasars are often found to be lensed. 
Furthermore, HzRGs are known to be associated with gi- 
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ant ellipticals in the local Universe (McLure & Dunlop 
2000). 

In this paper we present high-resolution observations 
of the CO J = 1 - emission from 4C 60.07 at z = 3.788 
using the Very Large Array^. Throughout we have as- 
sumed Ho = 65kms~^ Mpc~^, flM = 0.3 and = 0.7. 
In this cosmology the luminosity distance of 4C 60.07 is 
36.2 Gpc and 1" corresponds to 7.7kpc. 

2. 4C 60.07 

4C 60.07 is an ultra-steep-spectrum (USS) radio galaxy at 
redshift of z = 3.788 (Chambers et al. 1996; Rottgering 
et al. 1997). It displays a Fanaroff-Riley II (FR II) edge- 
brightened double-radio morphology (Fanaroff & Riley 
1974). The radio morphology of 4C 60.07 is seen in Figure 
n] which shows a VLA C-band (6 cm) archive image of 
4C 60.07. The system consists of two main bright hot 
spots separated by about 9". The south-eastern compo- 
nent is further comprised of two components (B and C). 
Continuum emission is also seen from the radio core which 
is located ^ 2" west- northwest of the C-component. 

The CO J = 4 — 3 line and 1.25 mm continuum 
emission from 4C 60.07 have been imaged using the 
IRAM Plateau de Bure Interferometer (Papadopoulos et 
al. 2000). The CO emission was found to emerge from 
two kinematically distinct components separated in ve- 
locity space by ~ 700kms~^. The component with the 
narrowest line profile (FWHM 150kms~^) is offset 
by ~ — 224kms^^ from the systemic velocity cor- 
responding to z — 3.788 and spatially coincident with 
the position of the radio core. The broader component 
{FWHM > 550 km s"^) peaks - 7" (- 30kpc) west of 
the radio core, and is offset by ~ 483kms~^ from the sys- 
temic velocity. Such large offsets between the redshift of 
the optical emission lines and the CO emission has been 
observed in several other high redshift systems, and is 
commonly attributed to the optical lines originating in 
strong outflows and winds (Guilloteau et al. 1999; Cox et 
al. 2002). Strong gravitational lensing can be ruled out as 
the origin of the double source and large apparent luminos- 
ity since not only do the two components have different 
line widths but they arc also offset in velocity with re- 
spect to each other. Furthermore, there are no indications 
from observations in the optical and radio suggesting that 
4C 60.07 might be lensed. 

In order to get a handle on the non-thermal contri- 
bution to the submm and CO fluxes, we used C and X 
band images from the VLA Archive to measure the ra- 
dio fluxes of the various components in 4C 60.07 at 6 cm 
and 3.6 cm. In Figure [21 we have plotted the submm/far- 
IR spectral energy distribution (SED) of 4C 60.07 along 
with the radio spectra of the A-component and the radio 

^ The Very Large Array (VLA) is operated by the National 
Radio Observatory, which is a facility of the National 
Science Foundation, operated under cooperative agreement by 
Associated Universities, Inc. 
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Fig. 1. 6 cm VLA map of 4C 60.07 overlaid on an I- 
band image obtained with the William Herschel Telescope. 
The radio hot spots are denoted by A, B and C. The 
cross marks the position of the radio core: R.A.(J2000): 
05''12™55?147, DEC.(J2000): +60°30'51'.'0. The insert in 
the bottom right corner shows the synthesized beam. 

core. The A-component (the dotted line in Figure|21) is the 
brightest component in the radio and furthermore has the 
shallowest spectral slope (ag cm™ — —1.4) which means 
it provides a strict upper limit on the non-thermal flux 
at submm wavelengths. The 1.25 mm continuum emission 
from 4C 60.07 is likely to be thermal in origin, i.e. from 
warm dust, since non-thermal processes are unable to ac- 
count for the observed flux, see Figure 2. Furthermore, the 
1.25 mm continuum emission is offset by ~ 4" to the west 
from the radio core position, and does not appear to be 
associated with the non-thermal emission. 

The radio core has a spectral slope of ag^^™ = —1.7. 
Extrapolating to the frequency of the CO J = 1 — line 

24 GHz) we find that the radio continuum is expected 
to contribute a non-neglible iS'24ghz = 0.03 mJy. In com- 
parison, the contribution from the thermal dust spectrum 
seems completely neglible, see Figure El 

3. Observations 

The CO J = 1 -Oline {lyrest = 115.27 GHz) from 4C 60.07 
is redshifted into the VLA's K-band (1.3-cm) receivers. 
Since non-thermal continuum emission from the radio core 
was expected, we used two IF pairs (right- and left-hand 
circular polarisation) centred on the two kinematically dis- 
tinct emission-line components, each pair set up such that 
one IF pair was centred on the line and the other pair 
was offset from this by 100 MHz to measure the contin- 
uum. The IF set-up for both the broad and narrow com- 
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Fig. 2. The submm/far-IR spectral energy distribution 
(SED) of 4C 60.07 (red line) along with the radio spec- 
tra of the A-component (blue dotted line) and the radio 
core (blue solid line). The submm points are taken from 
SCUBA observations by Archibald et al. (2001) while the 
IRAM 30m Telescope data-points at 252 GHz (1.25 mm) 
and 100 GHz (3 mm) are from Papadopoulos et al. (2000). 
The radio points are measured from C and X band data 
of 4C 60.07 from the VLA Archive. The positions and 
strengths of the CO J = 1 — and J = 4 — 3 lines rela- 
tive to the SED are shown in yellow (based on this paper 
and Papadopoulos et al. (2000)). 



ponent is detailed in Figure 13 The broad-line component 
was observed in continuum mode with a 50 MHz- wide IF 
centred at 24.035 GHz, which corresponds to a velocity 
coverage of 624 km s^^. Not only does this not properly 
cover the broad component but, since the observations 
were done in continuum mode, no information on the line 
shape was available, and one therefore has to rely on the 
IRAM CO J = 4 — 3 observations to infer a line width for 
the broad component. The narrow component was ob- 
served in spectral-line mode with a 7-channel IF centered 
at 24.095 GHz, each channel being 3.125 MHz wide. This 
was to avoid under-resolving the narrow line in velocity 
space. However, the line datasets from the 2002 December 
20 and 2003 March 09 are centred at 24.089 GHz, corre- 
sponding to a shift of two channels, in order to get a bet- 
ter spectral coverage at the low-frequency side of the line. 
Thus nine channels (denoted channels 1,..,9 in Figure 
covered the line, corresponding to a velocity coverage of 
350kms-^ 4C60.07 was observed in the VLA's CD, D 
and C configurations (see Table The D-configuration 
data, however, turned out not to have the spatial resolu- 
tion required to properly separate the (B,C) component 
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Fig. 3. Velocity and frequency coverage of the IF set-up 
used in our VLA CO J = 1 — observations. The dashed 
line outlines the systemic velocity corresponding to z = 
3.788. In red is shown the velocity-channels of the IRAM 
PdBI observations of Papadopoulos et al. (2000); channels 
in which CO J = 4 — 3 emission was detected are marked 
in thick red. 

from the radio core, and was therefore discarded. The lack 
of D-array data is very unlikely to result in " filtering out" 
any extended CO emission since such emission would have 
to be extended over scales > 60" (^ 460 kpc). The fact 
that we do not find any change in total flux after includ- 
ing the D-array data verifies this point. 

In total, after calibration overheads, we obtained 6.3 hr 
of integration time on source for both the broad and nar- 
row components. Calibration and data reduction was done 
using standard recipes in the NRAO AXVS Cookbook. 
The amplitude was calibrated with the quasars 3C 48 and 
3C286 at the beginning and/or end of each transit. The 
phase drift was calibrated using a fast-switching technique 
in which we observed the nearby source 04494-1-63322 ev- 
ery few minutes. The data were obtained in good weather 
conditions, and the rms phase fluctuations after calibra- 
tion was less than 10°. For the spectral-line mode observa- 
tions of the narrow component, the bandpass of the system 
was calibrated using 0319-1-415. In all cases the bandpass 
calibration averaged over the 7 channels was better than 
92%. Data taken at different times and in different con- 
figurations were combined using DBCON, weighting each 
dataset with the its total gridded weight. 



4. Results &. Analysis 

The radio morphology of 4C 60.7 in the broad and narrow 
component IF set-ups is shown in Figure 01 Figure 0^ and 
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Fig. 4. Naturally weighted, untapered maps, showing the morphology of 4C 60.07 in the two IF set-ups. Top row: 
The spatial resolution is 1'.'36 x 1'.'19 with the major axis P.A. = 89°. The rms noise is 0.030mJybeam~^. Contours are 
shown at -2, 2, 3, 4, 5, 6, 7, 8, 9, 10, 20, 30, and 40 x 35/xJybeam~^ Bottom row: The spatial resolution (FWHM) 
is 1"32 X 1'.'17 with the major axis P.A. = 275°, and the rms noise level is 0.051 mJybeam"^. Contours are shown at 
-2, 2, 3, 4, 5, 6, 7, 8, 9, 10, 15, and 20 x 60/iJybeam^^. The cross marks the position of the radio core (R.A.(J2000): 
05:12:55.147, DEC.(J2000): +60:30:51.0. 



Table 1. VLA Observations. 



Date 


Configuration 


tint /firs 


2001 Oct 2 


CD 


1.4 


2001 Oct 9 


CD 


1.7 


2001 Oct 11 


D 


1.7 


2001 Oct 19 


D 


2.8 


2002 Nov 12 


C 


1.7 


2002 Dec 20 


C 


1.5 



b show pure continuum emission and continuum emission 
plus the broad CO line emission, respectively; similarly 
for Figure^ and d but for the narrow component. 



In Figure which shows the continuum emission 
at 23.935 GHz, we recognise the two radio lobes seen in 
Figure 01 although we have failed to separate the B and C 
components from the radio core, let alone to resolve it into 
its two sub-components. The emission from the radio core 
is clearly stronger and more extended in IF2 which, in ad- 
dition to the continuum, contains the CO-emission from 
the broad component (Figure^) . A similar picture is seen 
for the narrow component in Figure^ and d. Here, the ra- 
dio core and B-component are separated in both IFs. The 
(large) negative spectral index (a ~ —1-7) of the radio 
core rules out the possibility that the excess emission seen 
in the IF2-maps is due to an increase in continuum emis- 
sion, since the IF2 maps are at a higher frequency than 
the IFl maps. The B and C radio hot spots have even 
larger negative spectral indices (Carilli et al. 1997) and 
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Fig. 5. a) Naturally weighted CO J = 1 — map of the broad component obtained by combining the entire data-set 
and tapering it with a Gaussian with a half-width at 0.30 amplitude of 200kA. The spatial resolution is 1'.'72 x 1'.'45 with 
the major axis P. A. =283°. The contours are at -2, 2, 3, 4, and 5 x cr, where a — 35 fi3yhea,m~^ . b) The CO J = 1 — 
emission tapered down to 60kA {FWHM — 3'.'33 x 3'.'12) overlaid as contours on a gray-scale representation of the 
J = 4 — 3 emission by Papadopoulos et al. (2000). The contours are -2, 2, 3, 4, 5, and 6(t with a = 0.4mJy beam^^, 
and the gray-scale range is 0.6 — 2.0mJy beam^^. 



can therefore not be the cause of the increased emission 
either. This suggests that what we are seeing in the IF2 
maps is continuum emission plus CO J = 1 — emission 
from the emission features detected in the CO J = 4 — 3 
maps. 

4.1. The Broad Line Component 

In order to disentangle the broad CO emission from the 
continuum, a box was put around the entire system and 
a CLEAN-component model of the continuum emission 
shown in Figure 0Jl was constructed. The CLEANing was 
stopped when the rms of the residuals reached the noise 
level in the image. This was then subtracted in uv-space 
from the IF2 data using the MPS task UVSUB, result- 
ing in uv-data from which maps free of continuum emis- 
sion could be produced. Figure shows the resulting 
CO J = 1 — map of the broad component. The lack of 
residual emission at the positions of the two radio lobes 
demonstrates the effectiveness of the continuum emission 
subtraction. The CO J = 1 — emission is detected at 5ct, 
peaking at a position which is coincident with the radio 
core position but with emission appearing to extend 2 — 3" 
to the west. In comparison, the broad component of the 
CO J = 4 — 3 emission was found to peak 4" south-west 
of the radio core (Papadopoulos et al. 2000). 

In Figure Eh we have overlaid the CO J = 1 — con- 
tours on top of a gray-scale image of the J = 4 — 3 emis- 
sion. The 1 — map has been tapered down to 60 kA in 
order to better match the resolution of the 4 — 3 data 
{FWHM = 8'.'9 X 5'.'5). The offset between the centroids 
of the 1 — and 4 — 3 emission is within the positional 
errors given the large J = 4 — 3 synthesized beam, and in 



general there is good spatial correspondence between the 
two. 

Assuming that the IF covers the entire line, we can es- 
timate the velocity-integrated J = 1 — flux density using 

/ Si^^^^dv = Ai^jF-^Sco, (1) 

J Ad '^obs 

where Avip is the width of the IF and Sco is the flux 
density. We find Sco = (0.27 ± 0.05) mJy which yields a 
velocity-integrated flux density of (0.15 ± 0.03) Jy kms~^, 
where we have used i/obs = 24.035 GHz and an IF band- 
width of 45 MHz (560kms-i) instead of 50 MHz due to 
bandpass rollover. The CO J = 1 — emission from the 
broad component in the uvtaper 200 map appears to be 
somewhat extended in the north-west direction, although 
the low signal-to-noise of the data does not warrant a firm 
conclusion on this issue. 

4.2. The Narrow Line Component 

The continuum subtraction in the case of the narrow- 
line emission was done in a similar fashion as for the 
broad component: a continuum model was constructed 
by combining all IFl channels and then subtracted from 
the IF2 channels. Since channels 1 and 2 had a differ- 
ent uv-coverage than channels 3 to 7 which again had 
a different uv-coverage than channels 8 and 9, three dif- 
ferent continuum models had to be constructed, one for 
each of these three sets of channels. In order to increase 
the signal-to-noise ratio we averaged neighbouring chan- 
nels during the imaging of the continuum-subtracted IF2 
channels. We then searched for any residual emission at 
the expected position of the narrow component and used 



6 



T.R. Greve et al.: Detection of CO J = 1 - in the z = 3.79 Radio Galaxy 4C 60.07 



1 .2 1 




-300 -200 

Velocity [km s '] 

Fig. 6. The spectrum of the narrow CO J = 1 — compo- 
nent in 4C 60.07. AU channels maps were tapered with 
a 200kA Gaussian. The rms noise in channels 1-2 is 
0.12mJy beam~^ at a spatial resolution of 1'.'4 x 1'.'3; the 
noise in channels 3-6 is 0.086 mJybeam"^ at a spatial 
resolution of 1'.'7 x 1'.'4, and the noise in channels 7-8 
is 0.09mJybeam~^ at a spatial resolution of 1'.'9 x 1'.'5. 
Above the line-profile we have plotted the velocity cover- 
age of the channels in which J = 4 — 3 emission was de- 
tected (Papadopoulos et al. 2000). The vertical dashed 
line corresponds to the systemic velocity (z = 3.788). 



the task imean to measure the flux. The resulting spec- 
trum of the CO J = 1 — emission is shown in Figure 
where we have also plotted the velocity-coverage of the two 
IRAM PdBI channels in which J = 4 — 3 emission from 
the narrow component was detected. Note that channels 
1-2 and 8-9 are somewhat more noisy than the other chan- 
nels, since the corresponding integration time is less. We 
detect CO J = 1 — emission in channels 3 to 6, which in 
velocity-space overlaps with the two IRAM PdBI chan- 
nels in which the strongest CO J = 4 — 3 emission is de- 
tected. Only very weak CO J = 4 — 3 emission was found 
in the thhd IRAM PdBI channel at 96.3095 GHz which 
is in agreement with the J = 1 — line profile. From a 
Gaussian fit to the line profile, we find the 1 — spectrum 
to peak at ~ (—220 ± 40)kms^^ offset from the sys- 
temic velocity of 4C 60.07 which is consistent with that 
found for the 4 — 3 line given the large velocity bins. The 
formal linewidth is AVfwhm — (165 ± 24) kms~^, again 
in good agreement with that of the CO J = 4 — 3 line pro- 
file. The fit yields a velocity-integrated flux density of 
0.09 ± 0.01 Jykms-\ see also Table 2. 

Combining channels 3 to 6 we obtain the velocity- 
integrated CO J = 1 — emission map shown in Figure 
7a. In Figure 7b we have overlaid contours of the 



CO J = 1 — emission on a gray-scale image of the 
J = 4 — 3 emission where the latter has been tapered 
down to 60 kA. The J = 1 — emission peaks at ^ 0'.'7 
east of the radio core which is consistent with the posi- 
tion of the narrow J = 4 — 3 component. Similar to what 
was found for the J = 4 — 3 line, in J = 1 — the narrow 
component appears to be more compact and less extended 
than the broad component. 

5. Discussion 

5.1. CO Luminosity and Molecular Gas Mass 

The observed CO J = 1 — line fluxes for the broad and 
narrow components in 4C 60.07 imply intrinsic CO lumi- 



nosities of L^q(1 — 0) 



(1.0 ± 0.2) X 10"Kkms-ipc2 
and (6.0 ± 0.7) x 10^" K kms^^ pc^, respectively. For the 
4-3 hne Papadopoulos et al. (2000) found Lco(4 - 3) = 
(7 ± 1) X 10i°Kkms-ipc2 and L^o(4-3) = (3.5 ± 
0.8) X 10^" Kkms^^ pc^ for the broad and narrow com- 
ponent, respectively, where we have computed the lu- 
minosities in the cosmology adopted here. From the 
CO J = 1 — line the molecular gas mass can be found 
using the well-known relation M(H2) = XcoLco(-'^ ~ 0) 
which relates the CO J = 1 — luminosity with the molec- 
ular gas mass (e.g. Strong et al. 1988). Xqo is the 
CO J = 1 — line luminosity to Il2-mass conversion fac- 
tor which in the extreme UV-intense environments found 
in local Ultra Luminous Infra-Red Galaxies, and pre- 
sumably also in high redshift galaxies such as 4C 60.07, 
has a value of about 0.8 (Kkms^^ pc^)~^M0 (Downes & 
Solomon 1998). In doing so we find molecular gas masses 
of M(H2) - 8 X 10^° Mq and M(F2) ~ 5 x 10^° Mq for 
the broad and narrow CO emitting components, respec- 
tively. Hence, even for a conservative, non-Galactic value 
of Xco we find that about ~ 10^^ Mq of molecular gas 
is associated with 4C 60.07. The estimated gas masses are 
in very good agreement with those of Papadopoulos et 
al. (2000). 

The total gas mass in 4C 60.07 will of course be 
larger once the neutral hydrogen has been accounted for. 
Assuming a value of M(HI)/M(H2) = 2 which is typically 
found in IRAS galaxies (Andreani, Casoli, & Gerin 1995) 
we find a total gas mass of Mgas — M(H2) + 2Af(H2) = 
2.4 X 10" Mq for the broad component and 1.5 x 10" Mq 
for the narrow component. In case that metal-poor gas is 
also present the gas mass can be even higher. 



5.2. FIR Luminosity and star-formation efficiency 

In the case where the far-IR luminosity is powered by 
a starburst and not an AGN, and all the stellar radi- 
ation is absorbed by dust, the far-IR to CO luminosity 
ratio, Lfir/ Leo, provides a rough measure of the inte- 
grated luminosity of massive stars responsible for heat- 
ing the dust {Lpiji) relative to the amount of fuel avail- 
able for star formation (LJ^q). We estimate the far-IR lu- 
minosity of 4C 60.07 to be Lpm ~ 4 x 10^^ , where 
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Fig. 7. a) The narrow component tapered at 200kA. The resolution is 1'.'72 x 1'.'42 at P.A. = 285°. Contours arc at -2, 2, 
3, and 4 x cr, where a = 55 /xJy beam"!, b) The CO J = 1 - emission tapered down to 60kA {FWHM = 3'.'22 x 3'.'07) 
overlaid as contours on a gray-scale representation of the J = 4 — 3 detection by Papadopoulos et al. (2000). The 
contours are -2, 2, 3, and 4cr with a = 0.6mJybeam~^, and the gray-scale range is 0.6 — 2.0 mJybeam"^. 



Table 2. Observed and Derived Properties for 4C 60.07. 



Parameter 4C 60.07 4C 60.07 

(broad) (narrow) 

R.A.(J2000) 05''12"'54f75 05''12'"55f30 

Decl.(J2000) +60''30'50'.'92 +60°30'52'.'29 

Scoci-o) [mJy] (0.27 ±0.05) (0.30 ±0.10) 

Scoci-ojAy [Jykms^^] (0.15 ± 0.03) (0.09 ± 0.01) 

Lco(i-o) [Kkms-^c^ (1.0±0.2) x lO" (6.0±0.7)xl0' 

M(H2) [Mq] 8 X 10^" 5 X 10^° 



we have used the 850- /im flux measurement of Archibald 
et al. (2001) and adopted a dust temperature and spec- 
tral index of Td = 50 K and /3 = 1.5, respectively. 
This yields a Lfir/Lco ratio of ~ (4 x IO^^Lq )/(1.6 x 
10"Kkms-ipc2) = 250Lo(Kkms-ipc2)-i, which is 
similar to values found in ULIRGs (e.g. Solomon et 
al. 1997). Carilli et al. (2002a) foimd eontinuum-to-line ra- 
tios of 350 and 323 in the QSOs BRI 1202-0725 {z = 4.70) 
and BRI 1335-0417 (z = 4.41), respectively. A somewhat 
lower value of 200 was foimd in the z = 4.12 QSO 
PSS J2322+1944 (Carilh et al. 2002b). 

The current star-formation rate measured in Mq yr~^ 
is given by 

SFR ~ {LfirW-^VLq )6mfSsb, (2) 

where Smf ~ 0.8 — 2 depends on the mass func;tion of the 
stellar population, and 633 is the fraction of the far-IR 
luminosity which is powered by the starburst and not the 
AGN (Omont et al. 2001). If we adopt 5mf = 0.8 and con- 
servatively assume that only 50% of Lfir is powered by 
the starburst {Ssb = 0.5) we find SFR ~ 1600 M© yr^^. 



Assuming that the total amount of molecular gas observed 
towards 4C 60.07 ends up as fuel for the starburst, we find 
that there is enough molecular gas to sustain the inferred 
star-formation rate for ~ 8 x lO'^yr. While this is short 
compared to the Hubble time it is still sufficient to produce 
a giant elliptical with a stellar mass of ~ 10"'^^ Mq . 

The efficiency with which stars are being formed, 
i.e. the rate of star-formation per solar mass of molec- 
ular hydrogen, is given by SFE = SFR/M(H2) or 
equivaiently LF/fi/M(H2). For 4C 60.07 we find SFE = 
LFiR/M{i{2) — 3001/0 Mq which is somewhat higher 
than the ~ 190^0 Mq ~^ reported by Papadopoulos et 
al. (2000) who derived their value based on their flux 
density measurement of Ssso/im = llmJy. Here we have 
used the measurement by Archibald et al. (2001) which 
yields a somewhat higher 850 /xm flux density of 17mJy, 
although still lower than the 22 mJy reported by Stevens 
et al. (2003). 

Local ULIRGs exhibit star-formation efficiencies com- 
parable to that of 4C 60.07, once the same Xco-factor has 
been used (Solomon et al. 1997). It is worth stressing that 
our detection of CO J = 1 — enables us to make a direct 
comparison of SFEs with that of local ULIRGs since the 
same gas mass measure (the J = 1 — line) is used for 
ULIRGs. 

The apparently high star-formation efficiencies found 
for the above systems, could be severely overestimated if 
the AGN contributes significantly to the far-IR luminos- 
ity. However, the detection of CO together with the fact 
that 4C 60.07 and other HzRGs appear extended on sev- 
eral tens of kilo-pc scales at submm-wavclengths (Ivison 
et al. 2000; Stevens et al. 2003) strongly suggests that the 
far-IR emission is powered by large-scale starburst and not 
the AGN. Here we must mention that while CO J = 1 — 
may be a good indicator of the total metal-rich H2 gas 
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resevoir, it may be a poor one regarding the dense gas 
phase {n > 10^ cm^'^) that "fuels" star-formation. The 
latter could be particular true in the tidally disrupted gi- 
ant molecular clouds (GMCs) expected in ULIRGs where 
a diffuse phase may dominate the CO J = 1 — emission 
but has little to do with star formation. This could be 
the reason why the Lpm/ L'qq ratio is found to be such a 
strong function of LpiR, increasing for the merger sys- 
tems associated usually with large far-IR luminosities. 
Interestingly, recent work shows that the SFE of dense gas, 
parametrised by the Lpj^/ Lhcn{i-o) ratio (the HCN 
J = 1 — critical density is 2 x lO^cm"^), remains con- 
stant from GMCs all the way to ULIRG system (Gao & 
Solomon 2003). 

5.3. Excitation Conditions 

From the detections of CO in 4C 60.07 we can infer the 
CO (4 — 3)/(l — 0) velocity/area-averaged brightness tem- 
perature ratio using r^^ = Lco(4-3)/-'^CO(i-o)' global 
line ratio, i.e. the line ratio obtained by combining the 
flux from the two components, is r^^ = 0.7^q'2. The line 
ratios for the broad and narrow components are 0.7^q'2 
and 0.6lo;2, respectively. Hence, given the large uncer- 
tainties we find no significant difference in the excita- 
tion conditions between the two components. Thus in this 
case H2 mass estimates solely from CO J = 4 — 3 assum- 
ing full thermalisation and optical thickness of the lat- 
ter (i.e. ~ 1) would not result in too large errors. 
However, the sub-thermal excitation of such high-J CO 
lines remains a possibility even in starburst environments 
(e.g. Papadopoulos & Ivison 2002; Greve et al. 2003). 

In Figure IHl we plotted the velocity- integrated CO line 
flux densities (normalised to J = 4 — 3) for the broad and 
narrow components in 4C 60.07 as well as the normalised 
hue fluxes for PSS J2322-I-1944. We have used a standard 
single-component large velocity gradient (LVG) code to 
interpret the observed line ratio. In fact, given the large 
range compatible with our measurements we can only put 
a lower limit on the density. Indeed, for both components 
the upper values for r^^ are compatible with LTE and op- 
tically thick emission where the ratio is no longer sensitive 
to the density. 

A lower limit on the average gas density can be 
set by the lowest possible value of — 0.4 allowed 
by the observations. Adopting a typical CO abundance 
A = Xco/{dV/dr) = 10-5(kms-ipc-i)-i, a Tcmb = 
(1 -I- z) X 2.73 K = 13.07 K and a Tk^„ = 50 K, which 
is a typical dust temperature in starburst environments 
(Colbert et al. 1999; Hughes et al. 1997), we find ^(Ha) > 
2 X lO^cm"^. Adopting a higher A = 10~^ (unlikely in 
such kinematically violent, UV- aggressive environments) 
lowers the aforementioned limit to ^ 600 cm^'^. A lower 
assumed temperature of Tkm = 30K does not change 
these lower limits by much. For the most likely value of 
r43 = 0.7 we find ^(Hz) 3 x lO^cm'^ (Tfc = 50K, 
A = 10"5(kms-ipc-i)-i). 



T 1 1 1 1 1 1 1 1 1 1 1 T-T-1 1 r 

o 
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Fig. 8. Velocity-integrated CO line flux densities from 
4C 60.07 are shown as filled squares. The integrated line 
fluxes have been normalised to the J = 4 — 3 line. Also 
shown are the line fluxes from the z = 4.12 QSO 
PSS J2322-I-1944 (filled triangles - CariUi et al. 2002b; Cox 
et al. 2002), the z = 4.69 QSO BRI 1202-0725 (open cir- 
cles - Omont et al. 1996; Carilli et al. 2002a), the local 
starburst galaxy M82 (open triangles - Mao et al. 2000), 
and the integrated emission from within the solar radius 
of the the Milky Way (filled circles - Fixsen et al. 1999). 
The red dotted line shows the line flux increasing as fre- 
quency squared which is expected for optically thick con- 
ditions. The green lines show results from a LVG-model 
with Tk^n = 50 K, Xco/{dV/dr) = 10"^ (kms"! pc'^)'^, 
and n(H2) equal to 3.0 x lO'^cm^'^ (sohd hue) and 2.2 x 
lO'^cm"'^ (dashed line). 



Comparing with the z = 4.12 QSO PSS J2322-f^l944 
which has a (4-3)/(l-0) line ratio of 1.4 ± 0.6 (CariUi et 
al. 2002b; Cox et al. 2002) it appears that the excitation 
conditions in 4C 60.07 are less extreme. Hence, these re- 
sults may indicate that the molecular gas in 4C 60.07 is 
not as dense as that seen in PSS J2322-I-1944. Such a com- 
parison, however, is prone to the effects of strong lensing. 
While 4C 60.07 is unlikely to be a strongly lensed sys- 
tem because of its CO components with widely different 
line widths, it is not the case for PSS 2322-1-1944 which is 
known to be lensed (Carilli et al. 2003). If the distribu- 
tion of the 1 — and J = 4 — 3 emitting gas in the source 
plane differs, differential magnification of the two lines can 
result in erroneous intrinsic line ratios. 
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5.4. What is the Evolutionary Status of 4C 60.07? 

A comparison between the molecular gas mass and the 
dynamical mass allows for the determination of the evo- 
lutionary status of a galaxy, while a comparison of its 
dynamical mass with that of present-day spiral or ellip- 
tical can point toward its possible descendant. Typically, 
dynamical masses are calculated assuming the gas is dis- 
tributed in a disk in Keplerian rotation (e.g. Carilli et 
al. 2002a). In the case of 4C 60.07, however, we have suffi- 
cient spatial and kinematical information to conclude that 
the two gas components do not belong to such a structure. 
The detection of two kinematically distinct gas resevoirs 
in 4C 60.07, each with a large gas mass, suggest a major 
merger event. Hence, a more plausible scenario might be 
that the two clouds are part of spherical system in the 
process of collapsing. Assuming the system is virialised, 
one can apply the virial theorem to derive the following 
expression for the dynamical mass 




where g is a factor which described the influence of non- 
gravitational forces on the virial equilibrium, and a ranges 
from 0.55 to 2.4 with a typical value of 1.5 which we 
shall adopt here (see Bryant & Scoville (1996) for details). 
A possible caveat to this argument comes from the fact 
that if the assumption of a virialised system is not true, 
eq. 13 will tend to overestimate the true dynamical mass. 
Adopting a radius of the sphere corresponding to half the 
separation between the two components, i.e. R ^ 15 kpc 
(corresponding to 2"), and a FWHM width equal to the 
velocity offset between the line centers of the two gas re- 
sevoirs, i.e. AVfwhm ~ 700kms~^, leads to an enclosed 
dynamical mass of Mdyn ^ 1-1 x 10^^ , where we 
have assumed that non-gravitational forces are negligible, 
i.e. q << 1. Note that since Bryant & Scoville (1996) de- 
rived eq.|31by examining a large number of clouds flUing up 
a spherical region, we have in the above assumed that the 
velocity and spatial separation between the two CO com- 
ponents is a fair approximation for AVpwHM f^nd R in 
eq.|31 Alternatively, we can do the more simple calculation 
of finding the virial mass required to keep the two clouds 
bound by a spherical potential. In this case the velocity 
dispersion becomes AV — (1/2)700 kms~^ = 350kms~^, 
and we find Mytnai = G^^RAV^ ~ 4.3 x 10" M© , where 
we have assumed R — 15 kpc. Projection effects will only 
tend to make R and AVfwhm, and therefore also the 
above mass estimates, larger. Thus, 4C 60.07 is a massive 
system with a dynamical mass comparable to that of gi- 
ant elliptical galaxies seen in the present day Universe. 
The inferred ratio of the molecular-to-dynamic mass for 
the system as a whole is M{H2)/Mdyn ~ 0.12, and the 
total gas fraction is Mgas/Mdyn ~ 0.35. The high gas frac- 
tion, along with the large velocity dispersion of the system 
suggest that 4C 60.07 is a very young system in its early 
stages of formation. 



The fact that such a massive {Mdyn ^ 10"'^^ Mq ) sys- 
tem has assembled at z — 3.8, which corresponds to a 
time when the Universe was only ~ 10% of its present age, 
seems to favour the monolithic collapse model (e.g. Eggen, 
Lynden-BeU & Sandage 1962; Tinsley & Gunn 1976) of 
massive ellipticals over the hierarchical formation scenario 
(e.g. Baron & White 1987; Baugh et al. 1996; Kauffmann 
& Chariot 1998). Further in support of the monolithic 
collapse scenario is the large star-formation rate [SFR ^ 
1600 M© yr~^) we infer from the far-IR luminosity. If such 
a large star-formation rate can be sustained it is capable 
of producing a giant elliptical in a time scale compara- 
ble to the dynamical time. The dynamical time-scale for 
a system with a mass M within a radius R is 

where tdyn is measured in years (see p. 37 of Binney 
& Tremaine 1987). For a giant elliptical with a mass 
Meii ~ 10^^ M© and a radius R = 30 kpc, we find 
tdyn — 1-2 X 10^ yr. Hence, in order to form a giant ellipti- 
cal within a dynamical time scale, star-formation rates of 
the order Meii/tdyn — 8000 M© yr~^ are required, which 
is comparable to what we find in 4C 60.07. It is important 
to note that even if the true dynamical mass or the star- 
formation rate of this system are less than the above esti- 
mates, the fact that at x = 3.8 about 1.3 x 10^^ M© in H2 
gas mass alone has assembled within 15 kpc points towards 
a major galaxy forming merger rather than a gradual ac- 
cumulation of mass over time and redshift. This mass can 
only be revised upwards and typically reach ~ 4x 10"'^"'^ M© 
if HI is also accounted for, thereby making the case for a 
monolithic collapse scenario even stronger. 

4C 60.07 is not the only system with extended CO 
emission. Molecular gas distributed on tens of kpc scales, 
sometimes in separate gas components, have been ob- 
served in high redshifts quasars (e.g. Carilli et al. 2002a; 
Papadopoulos et al. 2001) and radio galaxies (De Breuck 
et al. 2003; De Breuck et al. in prep.). Thus, there is evi- 
dence to show that at least some luminous active galaxies 
in the early Universe, in addition to harbouring super- 
massive black holes in their centres, are associated with 
massive reservoirs of molecular gas distributed on tens of 
kpc scales which are in the process of merging. This ob- 
served coevality between large scale mergers of gaseous 
subsystems and the epoch of AGN-activity might provide 
clues to the origin of the tight relationship observed locally 
between the velocity dispersion of spheroids and the mass 
of their central black holes (Ferrarese & Merritt 2000). 
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